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Abstract 

The control of chaotic systems impUes inducing an unpredictable system to follow a desired 
trajectory using the smallest "force". In low-dimensional continuous systems, one method is 
that of reconstructing the tangent space, so that the control may be concentrated on the most 
expanding directions. This scheme is hard to follow for high-dimensional (extended) systems. This 
is particularly true for extended systems that exhibit stable chaos, that is, systems which are not 
chaotic in the usual sense, but are unpredictable for finite perturbations. Prototypical of this class 
are cellular automata, aka completely discrete dynamical systems. Although usual indicators of 
chaoticity such as the maximum Lyapunov exponent may be defined for such systems, we show 
that the usual approach may lead to counter intuitive results, and that it is possible to exploit the 
characteristics of the system in order to reduce the distance between two replicas with less control. 

1 Introduction 

Control theory is a set of techniques for making a dynamical system behave in a desired way exerting a 
minimum effort. In particular, this technique may be applied to chaotic systems in order to make them 
follow a desired periodic orbit [1] or to synchronize a "slave" replica with a "master" one [2] . Most of 
the literature about control theory deals with low-dimensional systems modeled by a few differential 
equations. In such systems, the number of expanding directions is small and they behave smoothly in 
tangent space. 

In this paper we want to introduce the problem of controlling extended, highly non-linear dynamical 
systems. There is a class of systems, termed stable chaotic [3] which are nor chaotic in the usual 
sense of the sensitive dependence with respect to infinitesimal perturbations, but are nonetheless 
unpredictable. In particular, we shall concentrate on cellular automata (CA), which are discrete, 
deterministic dynamical systems. CA are widely used to model many systems in various fields, from 
computer science to earth sciences, biology, physics, sociology, etc. They are usually defined on a 
graph or a regular lattice, but may easily be extended to include mobile agents. The modeling of a 
system using cellular automata is conceptually much simpler than those using partial derivatives, and 
the evolution of such a system is easily performed by a digital computer, without rounding errors. 
However, for such systems continuity and smoothness (differentiability) do not apply. It is therefore 
hard to extend the usual techniques used in control thory and to define quantities like Lyapunov 
exponents and chaotic trajectories. 

In the case of master-slave synchronization, the "minimal strength" needed to synchronize a system 
is related to its chaoticity, defined by the lagest Lyapunov exponent in low-dimensional systems. For 



extended systems, the correspondence between the minimal strength and Lyapunov exponents may 
break down [4]. It is still possible to define derivatives of discrete systems [5], which prove useful in 
synchronization investigations [6]. 

In synchronization experiments, the "force" is generally applied blindly, without any relation with 
the dynamics. The corresponding synchronization effect is analogous to a directed percolation phase 
transition. The two systems synchronize when their difference goes to zero. Their difference grows due 
to their "chaotic" dynamics, along the directions identified by the Jacobian matrix of the evolution 
rule. The synchronization "pressure" reduces the paths along which a difference can propagate. When 
this reduction overcomes the chaotic growth, the system synchronizes. 

In control problems, one wants to exploit the knowledge about a system. It is therefore analogous 
to a synchronization problem of two different systems with a "targeted" force, that tries to "kill" the 
growing directions of the difference as soon as possible. We show how the concept of Boolean derivative 
and that of Boolean Jacobian matrix can be used to achieve this goal. 

This technique may be used also as a measurement technique: one system may be the experimental 
one, and the replica may be simulated on a computer. In this case, the two systems are in general 
different. We apply the control technique in order to synchronize two different cellular automata rules, 
one of which may be stochastic (i.e., influenced by external noise). 

2 Definitions 

Let us start our presentation by considering two smooth, chaotic maps 

x' = fix), 

y' = (1 +pf{x), 

where p is the control "strength" , / and g are two maps, x is vahicd at the discrete time t, and x/ 
is valued at i + 1 (the same applies to y and y'). The x is the "master", the y the "slave". The 
separation between both maps is u = x — y and the goal of control is to keep h = \u\ below a 
certain threshold (which may be zero) , using the minimum strength p. The two maps / and g may be 
different (for instance, they may use different parameters) or the same, in this case (synchronization) 
the synchronized state x = y is absorbing. 

If the desired trajectory is a natural one for the slave, control is equivalent to master-slave syn- 
chronization. 

x' = fix), 

y' = il-p)fiy)+pfix), 
and 

u' ^il-p){fix)- fiy)). 

For smooth maps, near the synchronization threshold Pc, it is possible to expand y(t) around the 
unperturbed trajectory xit), 

u' = il-p){fix)-fiy)):^il-p)^u. 

By iterating this map, one obtains the relation between synchronization threshold Pc and Lyapunov 
exponent A, 

= 1 - exp(-A). (1) 

The synchronized state is absorbing, since if for some time xit) = yit), then the control can be relaxed 
and the trajectories stay synchronized. However, for chaotic systems, this state is unstable. 
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2.1 Extended systems 

Natural systems, however, are rarely low-dimensional. We can extend the previous analysis by consid- 
ering a lattice of coupled maps, that may be thought as a stroboscopic view of a continuous system: 



x,{t + l) = f{g{x,.i{t),x,it),x^+,{t))). 



(2) 



where i = 1,. . . ,N. The function g represents the spatial coupling, it can be diffusive (linear) or highly 
nonlinear. The function / is the individual map, and can lead, when uncoupled, either to fixed points, 
stable cycles or chaotic oscillations. A perturbation may amplify exponentially in time by the action 
of /, but only linearly in time through the coupling (propagation to neighboring sites) . 

The dynamical properties of an extended system axe generally analyzed by means of the Lyapunov 
spectrum. Using vector notation, Eq. (2) can be written as 

x{t + l)=F{x{t)). 

The components of the Jacobian matrix of F are 



Jij{x{t)) 



i,j = 1,. . . ,N. For an infinitesimal perturbation 

d{t + l) = J{x{t))S{t). 

For instance, the Jacobian matrix of one-dimensional systems with nearest neighbor couplings has zero 
values except on the three central diagonals as 
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The eigenvectors of the Jacobian define the instantaneous tangent space of a dynamical system. The 
eigenvalues a, of the time product H^o "^{'^i^)) °f Jacobian matrices over a trajectory define the 
Lyapunov spectrum Aq > -^i > ^2 . . . with Aj = log(aj)/T [7]. It is generally assumed that a system 
is chaotic if Aq > and stable if Aq < 0. 

The largest Lyapunov exponent Aq (LLE) does not capture all the chaotic characteristics of an 
extended system. In general, a weak diffusive coupling reduces the LLE (since diffusion limits the 
exponential expansion along tangent space). Therefore, for small couplings, the maximum of chaoticity 
corresponds to uncoupled maps, but this situation may correspond to the easiest synchronizability (see 
Section 3.1). 



2.2 Stable chaos and cellular automata 

The scenario of extended systems may be more complex as we discuss below. Consider the map defined 
below and shown in Figure 1-a) 



f{x; a) 



'{&xY/2 0<a;<l/6, 

l-|6(l/3-a;)|°/2 1/3 < a; < 1/2, 

|6(a;-2/3)|"/2 1/2 < x < 5/6, 

,l-(6(l-a;))"/2 5/6 < a; < 1, 



(3) 
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1/6 1/3 1/2 2/3 5/6 1 
X 

(a) 



(b) 



Figure 1: (a) The plot of the map of (3) for different values of the parameter a. For a > 1, this 
map exhibits two attracting superstable fixed points, (b) Time evolution (downward) of a lattice of 
coupled maps. Color code: white=0, black=l, color=intermediate values. One can observe a transient 
disordered evolution (transient chaos) followed by a cellular automata pattern. 
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Figure 2: Possible inputs (xq, xi, X2), sum (t), output (/) and first-order derivatives of totalistic 
function R = IOIOI2 = lOjiQ. Since derivatives do not depend on configurations, it is a linear rule. 

This map is obviously stable for a > 1, with two fixed points = and Xi — 1, with interleaved 
basins that act as a sort of "frustration" when coupled, 

x.(t+l) = /( "-^('^ + "'^^)+"'+^('^ ), (4) 

i = 1, . . . , N with periodic noundary conditions. The system continues to be stable, exhibiting transient 
chaos (see Figure 1-b). After a transient, the eigenvalues of the Jacobian matrix go to zero, and all 
Lyapunov exponents go to —00. After a transient, the evolution is that of cellular automatn rule 150 
(see Section 3) , which is unpredictable for finite perturbations greater that 1/6 [8] . A similar behavior 
can be found in other continuous systems without direct correspondence to cellular automata [3]. 
Unpredictable stable systems are interesting since the synchronized state is stable. 

3 Cellular automata 

Cellular automata (CA) are completely discrete systems, defined as in Eq. (2), where Xi and / can 
assume values in a discrete set. In particular, we shall limit our study to Boolean CA, for which the set 
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of discrete values is {0, 1}. Since the function / is discrete, it can be defined by means of a complete 
enumeration of output given all possible inputs (look-up table) . We shall denote by r the size of the 
neighborhood, i.e., the number of cells whose state constitutes an input for the function /. Eq. (2) 
corresponds to r = 3. The case in which the function / is symmetric with respect to all inputs defines 
totalistic CA, since in this case one can consider that the value of the function / depends only on 
the sum of the values of sites in the neighborhood. While generic CA with range r are defined by 2^ 
entries in the look-up table, totalistic CA are defined by r -t- 1 entries. By arranging the output values 
of the look-up table as Boolean digits, one can compactly represent a CA rule as an integer number 
R, as shown for instance in Figure 2. 

Cellular automata may exhibit a large variety of dynamical behaviors. The number of possible 
states of a lattice of L Boolean cells is finite, and equal to 2^. Since the dynamics is deterministic, 
only limit cycles attractors are possible. One can divide the possible scenarios according with the 
number of attractors, the distribution of their basins and their period. For instance with r — trivial 
rules like rule have only one attractor with a large basin and period equal to one. The identity rule 
(which is not totalistic) has a large number of attractors (2^), each one with one state and period 
1. The majority rule IIOOI2 = 12|io has a intermediate number of attractors with short periods (1). 
"Chaotic" rules like rule IOIOI2 = 10|io exhibit cycles with very long period of the order of the total 
number of configurations as in Figure 3. Since in this case the period scales as an exponential of the 
size of the system, the difference between a periodic and aperiodic trajectory is not relevant (statistical 
quantities take similar values). Moreover, a defect or damage typically spreads in the configuration. 

It is possible [5] to extend the concept of derivative to cellular automata The Boolean derivative of 
F is the Jacobian matrix with components 

dFi{x) 

Ji,j = = Fi{xo, ...,Xj®l,.. .,XN-l) ® F.i{xo, ...,Xj,.. .,XN-l) 

{1 Fi changes when xj changes, ^ ^ 

Fi does not change when xj changes, 

where © denotes the sum modulo two. 

Many "standard" results may be extended to Boolean derivatives, for instance the Taylor expansion 



One can apply the "linear development" to discrete damages, and define a discrete Jacobian matrix, 
Eq. (5). Similarly to continuous systems, it is possible to define the largest Lyapunov exponent [6] 
related to the synchronization threshold (see Section 3.1). In contrast to continuous dynamics, defects 
can self-annihilate, so that the actual development of damage is different from the linearized one and 
they coincide only in the limit of vanishing damage, as shown in Figure 4. 
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Figure 3: Typical space-time patterns of "chaotic" rules. 
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damage spreading evolution in tangent space 
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Figure 4: Difference between damage spreading and evolution in tangent space 
3.1 Synchronization of extended systems 

There are many ways of "pushing" together two extended replicas. One possibility is "uniform" pushing 

y,it+l) = {^~p)F,iyit))+pF^{x{t)), 

for which the analysis presented above applies, with pc = 1 — exp(— Aq). This control is however quite 
difficult to be implemented experimentally in an extended system. Uniform synchronization of chaotic 
maps gives results similar to low-dimensional systems: pc — I ~ exp(— Aq). 
Another possibility is that of "pinching" synchronization 



y^it + l) 



Fiivit)) with probability 1 — p, 
Fi{x{t)) with probability p. 



In pinching synchronization, one has the possibility of applying the synchronization "strength" to a 
suitably chosen subset of sites. Pinching synchronization depends on coupling: uncoupled chaotic 
maps synchronizes for pc = 0. In general pc is larger for larger couplings [9]. 

In synchronization problems, synchronization is applied "blindly". In control problems, the goal 
is that of exploiting available information in order to apply a smaller amount of control (or achieve a 
stronger synchronization). 

4 Control of CA 

We study here the spatial application of synchronization 

{x' = F{x), 

I y' = (1 - p) F{y) © p F{x), 

where is the Hadamard (component by component) product, and the effect of synchronization 
Pi (z {0, 1} may depend on the position i. Therefore, the difference u evolves as 

u' = {l-p)Q{f{x)-f{y)) (6) 

and in the limit of vanishing distance, 

m'~(1-p)0Jm, (7) 
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(a) (b) (c) 

Figure 5: Plots of the different types of control for: (a) r = 3, i? = 1 (linear rule); (h) r = 3, R ~ 6 
(nonlinear rule); (c) r = 6, i? = 30 (nonlinear rule). For nonlinear rules, control 2 is worse and control 
3 is better than blind one (control type 1). 

where the scalar product Ju is computer modulo two. The control parameter is the average synchro- 
nization effort k — i^iPi)/N . The efficacy of synchronization (order parameter) is the asymptotic 
distance h = (YlnUi) /N . 

It is possible in principle to find the absolute minimum of k by computing the effects of all possible 
choices oi Pi, given an initial configuration x^^ = a;(0). This constitutes a great computational load. 
Since we are interested in possible real-time applications, we impose that the choice of = I may 
only depend on local information: the neighborhood configuration and a t = 1 time window. 

4.1 Implementation of control 

We investigate the following cases concerning how control p is applied: 

1. Blindly with probability p — k (standard pinching synchronization). 

2. With a probability p proportional to the sum of the first-order derivatives. 

3. With a probability p inversely proportional to the sum of first-order derivatives. 

In order to keep the implementation simple, instead of fixing k and computing the probability p, we 
let p be a free parameter, and measure the actual fraction of synchronized sites k and the average 
asymptotic distance h. The previous schemes only require information about x. If information about 
y or about the damage distribution u is available, the cost k is reduced by a factor h, since in this 
case we can apply the rule only when it is needed. 

4.2 Results 

Preliminary simulation results are presented in Figure 5. As expected, for linear rules there is no 
influence of the type of control, since all configurations have the same number of derivatives. For 
nonlinear rules, the observed behavior is the opposite of what is expected for continuous systems. 
Control 2, that minimizes the distance h for vanishing number of damages according to Eq. (7), 
gives worse results than the blind control 1. Control 3, inversely proportional to the sum of first-order 
derivatives, gives better results than the blind control 1. This result holds also for larger neighborhoods 
(Figure 5-c), but not for all rules. 

This surprising effect may be due to the fact that defects self- annihilate, as shown in Figure 6. In 
other words, we can exploit the characteristics of cellular automata (and other stable chaotic systems) 
in order to achieve a better control by exploiting the local contraction of the evolution rule. 
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Figure 6: Time evolution of defects for different types of control. Here r = 3 and R = 6, all cases 
starting from the same configuration. The effective probability p has been chosen so to have the same 
average control k in the three cases. One can notice that clusters of defects for control 3 are less dense 
than that of control 1 and 2. 

5 Conclusions 

Spatially extended stable systems (namely cellular automata) may exhibit unpredictable behavior 
(finite-distance chaoticity). The pinching synchronization threshold is related to this chaoticity. On 
the other hand, Boolean derivatives and discrete Lyapunov exponents may be used to characterize this 
kind of chaos. Synchronization may also be exploited for control in experimental situations. In the 
control problem one aims at discovering a protocol that keeps the distance h below a certain threshold 
with the minimum "effort" , given some constraints. We have chosen to investigate the behavior of two 
control schemes based on the local number of non-zero first-order derivatives, taking as reference the 
"blind" pinching synchronization protocol. 

We have shown that, differently from usual chaotic systems, one can exploit self-annihilation of 
defects to obtain synchronization with a weaker control, corresponding to the case in which the control 
is inversely proportional to the number of non-zero derivatives. 
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